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The  Institute  for  Defense  Analyse*  produces  three  kinds 
of  publication  for  distribution,  entitled  Report,  Study,  and 
Research  Paper. 

A  Report  embodies  the  results  of  a  major  research  project 
undertaken  by  IDA  and  is  intended  to  be  an  authoritative  con¬ 
tribution  on  its  subject. 

A  Study  is  a  less  formal  document  and  less  comprehensive 
in  scope  than  a  Report.  It  may  be  the  result  of  a  smaller  and 
more  narrowly  defined  research  project  or  it  may  be  a  sup¬ 
porting  technical  paper  prepared  in  connection  with  a  major 
project. 

A  Research  Paper  represents  the  work  of  one  or  more 
named  authors  but  is  subject  to  review  comparable  to  that  for 
publication  in  a  professional  journal. 
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FOREWORD 


The  work  reported  in  this  Study  is  part  of  a  continuing  effort  in 
the  analyses  of  alternative  civil-defense  systems  by  the  Institute 
for  Defense  Analyses  under  Contract  No.  OCD-OS-63-134  (dated 
June  28,  1963)  with  the  Office  of  Civil  Defense,  Department  of 
the  Army.  The  studies  are  being  performed  in  the  Economic  and 
Political  Studies  Division  of  IDA. 

This  Study  sets  forth  a  damage -limiting  strategy  for  allocating 
blast  and  fallout  shelters  to  the  inhabitants  of  a  city  or  metro¬ 
politan  area.  For  a  given  weapon,  a  method  is  outlined  for  ensur¬ 
ing  that  fatalities  within  any  local  area  under  study  do  not  exceed 
a  prescribed  proportion  of  total  population  irrespective  of  the 
actual  ground  zero  within  that  area.  Alternative  means  of  achieving 
this  objective  are  considered  and  a  least -cost  method  developed. 

The  proposed  population  protection  method  was  initially  con¬ 
ceived  in  June  1964  and  was  tested  that  summer  using  desk  calcula¬ 
tors.  The  results  substantiated  the  feasibility  of  the  strategy 
as  well  as  its  cost/effectiveness  potential. 

I  would  like  to  acknowledge  with  appreciation  the  contributions 
of  other  members  of  the  EPSD  staff  who  so  ably  assisted  in  this 
effort:  Martha  Strayhorn,  who  assisted  in  conducting  the  pilot  test 
of  the  shelter  allocation  procedure  last  summer;  Jane-Ring  Crane, 
John  D.  Wells,  and  Jane  Gleason,  who  developed  the  initial  computer 
model,  wrote  the  program  specifications  provided  in  Appendix  B  and 
supervised  the  production  runs  required  for  this  study  report;  and 
Phyllis  Weiner,  who  performed  all  of  the  technical  secretarial 
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duties  involved  in  the  development  and  publication  of  this  Study. 
Finally,  I  would  like  to  thank  Samuel  E.  Eastman,  the  Civil  Defense 
Project  Leader,  for  his  support  and  encouragement  throughout  this 
research  effort. 


Grace  J.  Kelleher 
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SUMMARY 


A  damage-limiting  strategy  for  allocating  blast  and  fallout  shelter 
protection  is  proposed  in  this  Study.  The  features  which  combine 
to  make  this  strategy  unique  are  its  relatively  fine-grained  local 
orientation  and  its  ability  to  meet  a  survival  percentage  criterion 
irrespective  of  the  actual  ground  zero  within  the  area  considered. 

Population  distributions  vary  from  city  to  city  and  must  be 
considered  in  detail--not  as  averages  or  aggregates — in  order  to 
provide  realistic,  optimal  protection  for  any  one  city.  The 
strategy  proposed  here  tailors  shelter  postures  to  the  conditions 
and  needs  of  individual  cities  or  local  areas.  This  local  approach 
could  be  used  to  develop  a  national  shelter  program  evaluating  the 
needs  of  many  cities  by  serial  application  of  the  shelter  allocation 
model . 

Shelter  postures  produced  under  this  strategy  consider  all 
potential  ground  zeros  within  the  protected  area  as  part  of  the 
shelter  allocation  process;  thus,  fatalities  from  immediate  blast 
effects  and  fallout  are  limited  to  a  stipulated  level,  irrespective 
of  where  an  assumed  weapon  might  be  delivered  within  a  target  city. 

Costs  are  minimized  in  the  shelter  allocation  process  by 
following  three  specific  decision  rules  described  in  Section  II. 

The  results  of  pilot  studies  reported  in  Section  III  confirm 
that  shelter  postures  developed  under  this  damage -limiting  strategy 
can  make  it  possible  to  limit  fatalities  to  a  stipulated  level1 
regardless  of  where  the  weapon  may  be  delivered  within  a  target 


1.  Subject,  of  course,  to  any  errors  involved  in  current 
estimates  of  basic  phenomena  and  effects  of  nuclear  weapons. 
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city.  These  studies  also  established  that,  by  judiciously  assigning 
shelters  of  various  strengths  to  different  areas  of  a  city,  the 
level  of  fatalities  can  be  limited  to  that  resulting  when  the 
strongest  shelters  under  consideration  are  provided  universally 
(at  much  greater  cost)  and  the  weapon  detonated  at  the  enemy's  best 
aimpoint.  For  example,  a  mixed-shelter  posture  comprising  300-  and 
100-psi  blast  shelters  and  fallout  shelters  having  a  pro¬ 
tection  factor^"  (PF)  of  40  or  greater  can  limit  fatalities  to 
2.2  percent,  given  a  10-Mt  surface  burst  on  the  Houston  at-home 
population.  This  posture  costs  only  $548  million,  yet  guarantees 
the  same  minimuifi 1  survival  level  as  a  posture  providing  300-psi 
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protection  to  everyone  at  a  cost  of  $1,014  million.  Similarly, 
for  the  same  attack,  a  posture  of  100-psi  blast  shelters  and  fallout 
shelters  with  a  PF  of  at  least  40  can  limit  fatalities  to  5  percent 
at  a  cost  of  only  $393  million.  A  posture  providing  100-psi  protec¬ 
tion  to  everyone  would  guarantee  the  same  minimum  survival  level 
at  a  cost  of  $562  million. 

The  shelter  allocation  model  is  described  in  Section  II  and 
specifications  for  computer  adaptation  are  provided  in  Appendix  B. 
Conclusions  of  this  Study  and  recommendations  for  further  develop¬ 
ment  of  the  model  are  set  forth  in  Section  IV. 


1.  The  protection  factor  of  a  shelter  is  a  measure  of  its 
ability  to  attenuate  radiation.  It  is  the  ratio  of  the  dose  which 
would  be  received  outdoors,  without  any  protection,  to  that  received 
at  a  particular  location  in  a  structure. 
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INTRODUCTION 


Earlier  civil-defense  studies  have  been  concerned  primarily  with 
estimating  casualties  from  a  generalized  attack  upon  the  United 
States,  based  upon  varying  assumptions  of  population  protection  in 
different  geographical  areas.  These  national  studies  as  well  as 
the  relatively  few  local  studies  have  assumed  that  the  enemy 
would  target  optimally  and  that  delivery  accuracy  could  be  approxi¬ 
mated  by  random  error  estimates  in  the  study  models.  This  method 
has  usually  precluded  meaningful  comparison  of  different  study 
results  because  of  differences  in  the  actual  ground  zeros  employed. 
More  important,  the  results  of  such  studies  leave  unanswered 
a  number  of  rather  compelling  questions :  Suppose  the  weapon  is 
not  optimally  delivered.  Could  it  not  detonate  at  ary  point  in 
the  city?  Should  our  defenses  be  keyed  to  estimates,  or  perhaps 
mere  conjectures,  of  the  enemy’s  targeting  intentions,  the  reliabil¬ 
ity  of  his  weapons,  and  his  delivery  accuracy? 

The  problems  posed  by  ground  zero  uncertainty  are  focused  more 
clearly  by  considering  the  full  range  of  future  potential  enemies 
and  their  relative  achievements  in  the  development  and  testing  of 
nuclear  weapons  systems,  i.e.,  the  n-th  country  problem.  For 
example,  the  center  of  the  business  district  of  a  city  might  be  the 
best  target  for  an  attack  against  population.  A  shelter  posture 
could  be  designed  to  counter  effectively  a  weapon  delivered  near 
that  point.  But,  owing  to  retargeting  or  error,  the  weapon  could 
detonate  in  a  densely  populated  suburban  area  within  the  city  and 
kill  almost  as  many  people  as  if  there  were  no  shelter  posture  at 
all.  The  shelter  system  would  have  bought  very  little  protection, 
if  any. 
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The  shelter  strategy  proposed  in  this  Study  circumvents  the 
uncertainties  involved  in  weapons  delivery;  fatalities  from  imme¬ 
diate  blast  effects  and  fallout  are  limited  to  a  given  level,  no 
matter  where  the  weapon  were  detonated  in  the  target  area.  An 
obvious  way  to  ensure  this  outcome  would  be  to  allocate  the  highest 
level  of  protection  to  everyone  in  the  potential  target  area.  Such 
a  posture  would  limit  fatalities  to  the  level  expected  from 
a  weapon  delivered  at  optimum  ground  zero.  A  number  of  such  univer 
sal  postures  were  assessed  earlier  and  were  found  not  only  to  over- 
protect,  but  also  to  entail  prohibitive  costs.  Clearly,  a  more 
efficient  technique  was  needed  to  meet  damage -limiting  objectives. 

A  shelter-allocation  strategy  that  is  both  damage-limiting 
and  least-cost  has  been  developed  and  is  described  in  this  Study. 
The  features  which  combine  to  make  this  strategy  unique  are  its 
relatively  fine-grained  local  orientation  and  its  consideration  of 
all  potential  ground  zeros  in  the  development  of  defensive  shelter 
postures.  Because  of  this  latter  feature,  the  postures  produced 
by  this  strategy  can  make  it  possible  to  limit  or  upper-bound 
fatalities  to  a  stipulated  level,  independent  of  ground  zero. 

A  model  for  executing  the  strategy  will  be  presented  along  with 
specifications  for  computer  adaptation.  As  will  be  demonstrated, 
postures  developed  by  this  strategy  against  a  single  weapon  of 
a  given  yield  will,  at  the  same  cost,  also  prove  highly  (if  not 
optimally)  effective  against  weapons  of  different  yields. 
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II 


SHELTER  ALLOCATION  MODEL 

The  proposed  strategy  is  to  provide  a  given  city  or  metropolitan 
area  the  least-cost  shelter  posture  that  will  limit  fatalities  from 
a  given  weapon  to  a  stipulated  upper  level  (a)1  regardless  of  where 
the  weapon  is  delivered  in  the  area.  The  strategy  is  executed  with 
a  unique  shelter  allocation  model.  The  narrative  description  of 
the  model  below  (IIA  and  B)  is  followed  by  a  detailed  discussion 
of  the  inputs  and  assumptions  (IIC)  and  the  computation  procedure 
(IID) . 

A.  DESCRIPTION  OF  THE  MODEL 

In  preparing  for  shelter  allocation,  the  model  determines  the 

distances  from  ground  zero  at  which  shelters  of  each  strength  used 

in  the  posture  would  be  required  for  protection.  The  regions  of 
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interest  about  ground  zero  for  a  single  10  Mt- surface  burst,  give 

a  choice  of  blast  shelters  (100-psi  or  35-psi)  and  fallout  shelters 

(PF  =  40  and  a  7.3-psi  blast  rating),  are  shown  in  Figure  1,  which 

2 

is  a  matrix  of  1-km  elements.  Having  established  these  regions  of 
effectiveness,  the  model  evaluates  all  potential  ground  zeros^  and 
allocates  shelters  as  required  to  limit  fatalities  at  the  stipulated 
upper  level,  a. 


1.  A  glossary  of  notation  is  included  as  Appendix  A. 

2.  A  0.50  fission  ratio  was  used  in  calculating  the  effects  of 
all  weapons  discussed  in  this  Study. 

3.  Each  element  of  the  population  matrix  is  considered  to  be  a 
potential  ground  zero.  In  effect,  the  template  depicted  in  Figure  1 
is  repeatedly  shifted,  so  that  the  ground-zero  square  falls  in  turn 
on  each  kilometer- square  element  of  the  population  matrix,  and  for 
each  new  ground  zero  the  number  of  blast  fatalities  is  calculated. 
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Legend: 


Region  Overpressure  (o)  Adequote  Blast/Follout  Protection 

HHmHj  =  Rq  o  >  100  psi  None 

llillllll  =  R|  35  <  o  s  100  psi  100  psi  shelter 

I  I  I  I  -  R2  7.3<  o  s  35  psi  35  psi  shelter 

j  |  1  |  =  R3  o  £  7.3  psi  Fallout  shelter,  PF  =  40 

Eil  =  Ground  Zero  ScalerQ  =  1  km3 

FIGURE  1  A  Basis  for  Determining  the  Protection  Required  to  Survive  at  Critical 

Distances  from  Ground  Zero,  10  MT  Surface  Burst 

(Given  a  choice  of  100  or  35  psi  Blast  Shelters  or  fall- 
,'M'1  out  protection  in  light  steel  frame  bldgs.) 


In  evaluating  each  ground  zero,  the  model  first  determines  how 
many  people  are  within  the  sure-kill  area  (R^  in  Figure  1).  This 
number  is  compared  with  the  stipulated  upper  limit  on  fatalities  to 
determine  the  number  of  fatalities  allowable  in  other  regions  about 
ground  zero.  For  example,  if  the  fatality  level  stipulated  for  the 
posture  is  50,000  persons  and  30,000  were  located  within  the  sure-kill 
region,  only  20,000  fatalities  are  allowable  elsewhere.  That  is,  all 
except  20,000  of  the  people  outside  the  sure-kill  radius  must  be 
adequately  protected  against  blast  and  fallout  effects.  The  next 
problem-  is  to  decide  which  .20,000  persons  will  not  be  protected. 

Sine  a  the  objective  is  to  minimize  the  cost  of  the  required  posture, 
inis  decision  is  based  on  the  costs  per  person  protected.  The  matrix 
elements  with  the  highest  protection  costs  per  person  are  identified 
and  left  unprotected,  while  the  rest  of  the  elements  are  allocated 
the  shelter  required  in  their  particular  region  about  ground 
zero.  In  evaluating  succeeding  ground  zeros,  all  shelters  previously 
allocated  are  applied  as  assets,  and  additional  protection  is  pro¬ 
vided  only  as  necessary  to  hold  fatalities  to  the  stipulated  level. 

After  all  ground  zeros  have  been  evaluated  for  appropriate 
blast  protection,  fallout  shelters  are  allocated  to  those  elements 
where  blast  shelters  were  not  previously  assigned.  The  final 
posture  is  then  costed. 

B.  SPECIAL  FEATURES  OF  THE  MODEL. 

The  principal  features  which  combine  to  make  this  model  unique 
are  its  local  orientation  and  its  deterministic  consideration  of 
all  potential  ground  zeros  in  the  development  of  shelter  postures. 

It  also  applies  a  minimum-cost  approach  to  the  problem.  Each  of 
these  characteristics  is  discussed  separately  below. 

1.  Local  Orientation 

Shelter  postures  are  tailored  in  this  model  to  the  conditions 
and  needs  of  individual  cities  or  local  areas.  This  local  approach 
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could,  of  course,  be  used  to  develop  a  national  shelter  program 
evaluating  the  needs  of  many  cities  by  serial  application  of  the 
model. 

A  basic  input  to  the  shelter  allocation  model  is  a  target 
value  matrix  in  which  the  population  of  the  target  city  is  broken 
down  according  to  uniform  subdivisions;  one -kilometer- square  matrix 
elements  are  used  in  the  current  version  of  this  model.  The 
assumed  time  of  attack  determines  whether  the  at-work,  at-home  or 
rush-hour  population  in  each  grid  element  will  be  the  basis  for 
shelter  assignments. 

2.  Deterministic  Consideration  of  All  Potential  Ground  Zeros 

All  grid  elements  are  considered  to  be  potential  ground  zeros. 
Each  one  is  separately  evaluated  in  order  of  target  interest. 

Thus,  this  is  a  deterministic  as  opposed  to  a  probabilistic  or 
"risk"1  model.  Shelter  is  allocated  as  necessary  to  ensure  that 
fatalities  from  the  assumed  weapon  will  not  exceed  the  stipulated 
level,  regardless  of  where  the  weapon  is  delivered  on  the  target 
city. 

3.  Minimum-Cost  Approach 

Three  decision  rules  for  minimizing  costs  are  applied  in  the 
shelter  allocation  process. 

(1)  Shelters  of  a  given  strength  are  not  allocated 
unless  essential  to  maintaining  the  stipulated  survival 
level. 

(2)  When  additional  protection  is  required,  shelters 
of  a  given  strength  are  allocated  only  for  those 
elements  in  which  they  are  both  necessary  and  effective 
for  protection  against  the  specific  ground  zero  being 
evaluated. 


1.  In  a  ’’rnsk”  .v.odol,  probabilities  would  be  assigned  to 
a  number  of  potential  ground  zeros  and  the  average  or  expected 
damage  calculated  for  each  area  affected.  In  contrast,  the  pro¬ 
posed  model  deals  w?.th  the  maximum  fatalities  that  could  be 
incurred  from  a  given  weapon. 
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(3)  Cost  dictates  the  choice  of  specific  elements  to  be 
sheltered.  Priority  is  given  to  those  elements  that  can 
be  protected  at  the  lowest  average  cost  per  required 
shelter  space. 

The  element-by-element  differences  in  cost  per  required  space 
are  illustrated  in  Appendix  C. 

Once  shelter  of  a  given  strength  has  been  allocated  to  a  popu¬ 
lation  element,  it  may  not  be  downgraded.  However,  it  may  be  re¬ 
placed  with  a  stronger  shelter  if  the  process  of  limiting  fatalities 
for  another  potential  ground  zero  so  dictates. 

i  I 

The  shelter-allocation  model  provides  a  feasible  minimum-cost 
solution  to  the  problem.  Whether  it  produces  the  minimum  solution 
will  not  be  argued  at  this  stage  of  the  model's  development.  In 
seeking  such  a  solution,  linear  programming  techniques  were  con¬ 
sidered.  However,  their  application  was  necessarily  deferred 
because  of  the  prohibitive  number  of  constraints  now  associated 
with  the  problem.  Available  linear  programming  models,  e.g. , 

LP-90,  will  not  accommodate  the  total  constraints  involved.  Poten¬ 
tial  means  of  reducing  the  constraints  include: 

(1)  reducing  the  number  of  grid  elements  by  increasing 

their  size  (presently  one-kilometer  squares);  ^ 

(2)  reducing  the  total  area  of  the  target  value  matrix, 
thereby  reducing  the  total  number  of  elements  involved; 

(3)  using  fewer  sizes  and  strengths  of  shelter  (perhaps 
some  ideal  shelter  complex);  or 

(4)  considering  the  present  cost/effectiveness  approach 
as  the  first  step  toward  the  minimum-cost  solution  and 
then  using  the  results  to  determine  those  areas  of  the 
population  matrix  to  which  an  available  linear  program 
could  be  applied  in  an  attempt  to  obtain  a  cheaper 
solution.  The  areas  appear  to  be  those  where  a  concen¬ 
tration  of  one  level  of  protection  ends  and  another 
begins. 

The  principal  difference  between  a  linear  programming  procedure 
and  that  being  recommended  here  is  that  the  former  would  permit  simul¬ 
taneous  consideration  of  all  matrix  elements  as  potential  ground  zeros. 
In  the  present  procedure,  elements  are  evaluated  one  by  one,  in  a 


7 


prescribed  sequence,  crediting  as  assets  all  shelters  allocated 
earlier  in  the  sequence. 

C.  INPUTS,  GUIDELINES,  AND  ASSUMPTIONS1 

The  required  inputs  to  the  shelter  allocation  model  are 

(1)  target  value  matrix; 

(2)  the  threat; 

(3)  shelter  choices--sizes  and  strengths; 

(4)  shelter  costs; 

(5)  the  level,  a  ,  at  which  fatalities  are  to  be  upper- 
bounded;  andn 

(6)  the  psi  level  above  which  blast  shelters  would  be 
required  for  population  protection. 

Each  of  these  is  discussed  separately  below. 

A  target  value  (population)  matrix  is  input  for  the  target  city 
to  be  protected.  This  matrix  was  described  earlier  on  pages  5  and  6. 

The  threat  to  be  countered  by  the  shelter  posture  is  also  an 
input.  This  includes  the  yield  of  the  weapon  and  height  of  burst. 

Wind  direction  and  speed  are  considered  externally  in  determining 
the  level  of  fallout  protection  required. 

The  shelter  strengths  and  sizes  available  for  allocation  and  their 
respective  costs  are  also  inputs  to  the  model.  See  Appendix  C  for 
examples  and  cost  sources. 

Shelter  costs:  The  total  population,  p^.,  of  a  grid  element 
(identified  in  aggregate  only)  is  the  basic  decision  unit  affecting 
cost.  Each  p^..  will  be  provided  a  uniform  level  of  protection,  i.e., 
shelters  of  the  same  psi  rating,  or  the  same  fallout  protection 
factor,  for  everyone  in  the  element. 

The  level  (o^)  at  which  fatalities  are  to  be  upperbounded  is 
also  an  input.  A  series  of  postures  can  be  developed  and  costed 


1.  Appendix  A  is  a  glossary  of  all  notation  used  here  and  in 
Section  IID.  A  set  of  specifications  for  computer  adaptation  of  the 
model  is  provided  in  Appendix  B. 
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-'based  on  alternative  levels,  thus  providing  discrete  points  for 
the  construction  of  cost/effectiveness  curves.  The  range  of 
feasible  <*n  levels  is  determined  by  the  threat  to  be  countered  and 
the  alternative  shelter  strengths  available  for  the  posture. 

(1)  The  minimum  level  at  which  fatalities  can  be  upper-bounded 
(cm)  is  the  number  of  people  that  would  be  killed  by  blast 
effects  from  the  stipulated  weapon  if  the  strongest  shelter 
available  for  the  posture  C s^)  were  provided  universally. 

That  is,  even  if  s,  shelters  were  allocated  to  everyone,  the 
assumed  weapon  could  kill  as  many  as  a,  people  from  blast 
effects,  depending  upon  its  actual  ground  zero  within  the  city. 

(2)  The  maximum  level  of  fatalities  for  which  at  least  some 
Sn  shelter  would  be  required  (am)  is  one  less  fatality  than 
the  upper  limit  protected  by  the  next  strongest  shelter  (S2). 

For  example,  if  fatalities  could  be  limited  to  50,000  with 

s2  blast  shelters  (or  with  fallout  shelters  if  blast  shelters 
or  only  one  strength  are  being  allocated),  the  maximum  a 
requiring  allocation  of  any  s,  shelter  is  49,999.  If  any 
higher  a  were  stipulated,  the  model  would  not  allocate  any  s-, 
shelter,  as  it  would  not  be  required.  (The  model  will  not 
allocate  a  stronger  shelter  if  a  less  expensive  shelter  of 
lower  strength  will  provide  an  equivalent  number  of  survivors. 
Likewise,  it  will  not  allocate  any  blast  shelters  if  fallout 
shelters  will  ensure  the  stipulated  number  of  survivors.) 

The  fatality  levels  corresponding  to  and  need  not  be  specified 

in  advance,  as  their  values  will  be  computed  as  a  preliminary  step 

in  the  shelter  allocation  procedure  itself. 

The  minimum  psi  level  for  which  blast  shelters  would  be  required 
for  population  protection  is  also  an  input  and  can  be  varied. 

D.  COMPUTATION  PROCEDURE 

The  first  three  steps  in  the  computation  are  in  preparation  for 
the  shelter  allocations;  the  fcurth  and  fifth  are  the  allocation 
procedure,  and  the  sixth  is  the  costing  procedure. 

Step  1:  Establish  Shelter  Effectiveness  Radii 

In  preparing  for  shelter  allocation,  it  is  necessary  to  determine 
the  distances  from  ground  zero  at  which  a  shelter  of  specified  strength 
k  will  be  effective.  The  distance  to  be  calculated  is  r^,  the  outer 
limit  of  the  region  R^  within  which  a  shelter  of  a  given  strength  in 
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the  posture  will  protect  and  a  weaker  shelter  in  the  posture  will  not. 
Distance  r^  is  determined  from  overpressure -versus-distance  formulas,1 
properly  scaled  for  weapon  yield  and  height  of  burst.  Figure  1  depicts 
for  a  10-Mt  surface  burst,  the  region  of  effectiveness  for  each  of 
three  shelter  strengths:  100-  and  35-psi  blast  shelters,  and  fallout 
shelters. 

Step  2:  Determine  the  Minimum  Level  at  Which  Fatalities 
Can  Be  Upper  Bounded 

PQ  is  the  population  in  the  unprotectable  region  RQ,  with  outer 
radius  rQ,  about  a  given  ground  zero.  In  this  region  the  over¬ 
pressure  from  the  assumed  weapon  will  exceed  the  psi  rating  of  even 
the  strongest  shelter.  RQ  contains  all  matrix  elements  located  at 
distances,  d^. ,  from  ground  zero  that  are  less  than  or  equal  to  rQ. 

The  distance,  d..,  of  an  element  in  the  i-th  row,  i-th  column  from 
ij 

a  ground  zero  in  the  m-th  row,  n-th  column  is  determined  by  the 
following  formula: 

d^.  =  V  (i-m)^  +  (j-n)2 


If  d..  ^  r  ,  the  element  i,j  is  included  in  R  ,  and 
in  o  o 


Po=S 


*0 
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Calculate  P  for  each  grid  element,  considering  that  each  is  a 
°  2 

potential  ground  zero.  The  highest  PQ  calculated  is  a the  minimum 
level  at  which  fatalities  can  be  upper-bounded  with  the  shelter  choices 
at  hand. 


1.  Based  on  data  provided  in  Figure  3.66  of  S.  Glasstone  (Ed.), 
The  Effects  of  Nuclear  Weapons,  U.S.  Department  of  Defense  (U.S. 
Atomic  Energy  Commission:  Washington,  D.  C. ,  April  1962),  p.  135. 

2.  As  will  be  described  later,  the  resultant  P  values  also 
provide  the  basis  for  ordering  elements  for  evaluation  as  potential 
ground  zeros  in  the  shelter  allocation  process. 
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Step  3:  Determine  g  ,  the  Highest  Fatality  Level  that  Would  Require 
Any  Shelter  Allocation 

is  the  total  population  in  regions  RQ  and  about  a  given 
ground  zero. 


pi =  re  pij +  f  pij  • 

R0  K1 


The  second  strongest  shelter  under  consideration  (s2)  will  not 
protect  within  R^  (bounded  by  radius  r^),  but  will  protect  outside 
r^.  Therefore,  to  bound  fatalities  at  would  require  no  s^ 
shelters.  The  populations  of  all  elements  with  d^  <  r^  are 


included  in  P^. 

Calculate  P^  for  each  grid  element  considering  that  each  is 
a  potential  ground  zero.  The  highest  P^  calculated,  Pjjnax*  t*ie 
minimum  level  at  which  fatalities  can  be  upper  bounded  with  s2 
shelters.  Thus, 

Plmax  1  =  am  * 

the  highest  fatality  level  that  would  require  allocation  of  any 
s^  shelter. 

Step  4 :  Determine  Shelter  Requirements  and  Allocate  Shelter 

Array  the  blast-shelter  choices  in  order  of  psi  rating,  highest 
to  lowest.  Classify  the  strongest  shelter  s^,  the  next  strongest 

s2,  etc.  Allocate  shelters  one  strength  category  at  a  time,  in 

,  1 
the  order  (k  =  1,  2,  3,  ...).  To  allocate  s^  shelter: 


1.  An  alternative  for  this  step  ^s  discussed  at  the  end  of 
this  Section.  See  page  14. 
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j  Step  4a:  Array  the  matrix  elements  from  highest  to  lowest 

according  to  the  numbers  of  persons,  V.  vulnerable  even  with  s, 

mn  k 

shelters,  to  determine  the  order  in  which  they  are  to  be  evaluated 
as  potential  ground  zeros. ^ 

For  the  strongest  shelter,  s^, 

Vmri  =  P0  . 
mn  o 

The  PQ  values  were  calculated  in  Step  2  above. 

For  other  shelter  strengths,  s^  (k  =  2,  3,  ...), 

Vnm  =  Po  +  <U1  +  U2  +  •••  +  5 

where  (k  =  1,  2,  3,  ...)  is  the  total  population  in  those  ele¬ 
ments  of  R^  which  were  not  allocated  s^  shelters. 

Step  4b:  Evaluate  each  element  a^  as  a  potential  ground 
zero  and  allocate  shelter  as  follows : 

Calculate  f^,  the  allowable  fatalities  in  R^,  where 

fk  =  an  -  Vmn 

(an  is  an  input,  as  described  earlier  in  this  Section  on  pages 
7  and  8). 


Calculate  U,  where 


uv  =  r  u..  ; 


and  the  u^j  are  the  populations  in  elements  i,j  which  have  not 
already  received  s^  or  strongsr  shelter. 


1.  Ground  zeros  could  be  evaluated  inkany  conceivable  sequence 
and  still  accomplish  the  objectives  of  the  proposed  strategy.  Sim¬ 
ultaneous  consideration  of  all  ground  zeros,  using  a  linear  pro¬ 
gramming  approach,  might  provide  a  somewhat  lower-cost  solution  as 
discussed  in  IIB3. 
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If  <  f^>  no  additional  s^  shelter  is  required  for  this 
ground  zero. 

If  Uk>  f^>  some  s^  shelter  is  required.  Calculate  the 

average  cost  per  required  shelter  space,  c^^,  highest  to  lowest. 

This  array  will  be  used  to  determine  which  elements  are  to  be  left 

unsheltered  in  at  this  point  in  the  allocation  process. 

Starting  with  t*e  unprotected  element  having  the  highest  c^^ 

in  R^,  select  to  leave  unprotected  each  element  down  the  array 

until  either  all  u.  .’s  have  been  selected  and  the  sum  <  f,  or  the 

13  -  k  — 

sum  of  the  u.  .’s  equals  or  just  exceeds  fv.  If  this  latter  sum 

1 J  K 

just  exceeds  f,  ,  reduce  it  by  the  last  u..  added.  The  u..  elements 
remaining  in  the  sum  are  left  unsheltered  and  all  others  in  R^  are 
assigned  s^  shelters  in  the  numbers  and  sizes  required. 


Step  4c:  All  elements  other  than  the  first  are  evaluated  in 
sequence,  highest  to  lowest  V^.  Repeat  Step  4b  above. 

Step  4d;  Repeat  Steps  4a  through  4c  for  each  s^  shelter  until 
all  choices  (k  =  1,  2,  3,  ...)  have  been  considered  and  allocated 
as  necessary. 


Step  5 :  Allocate  Fallout  Shelters 

When  all  shelter  choices  s^  s2,  ...,  sn  have  been  allocated 

in  accordance  with  the  above  procedure,  the  blast  shelter  portion 

of  the  posture  is  completed.  The  next  step  is  to  allocate  adequate 

fallout  shelter  to  all  the  p.  .’s  which  were  not  assigned  special 

1 

blast  protection. 

Step  6 :  Costing 

Cost  this  posture  using  available  cost  factors  for  the  shelter 
choices  (strengths  and  sizes)  allocated.  The  rr^st  economical  com¬ 
bination  of  shelter  sizes  that  will  protect  a  given  p..  with  the 


1.  The  required  level  of  fallout  protection,  e.g.,  PF  >  40, 
for  the  given  threat,  is  determined  by  a  separate  IDA  computer 
routine  called  FALDIK. 
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assigned  shelter  strength,  s^,  is  always  used  as  the  basis  for 
blast  shelter  costing.  In  the  present  model,  required  fallout 
shelter  protection  is  costed  by  summing  the  unsheltered  population 
at  the  end  of  Step  3  above,  and  then  applying  a  single,  average 
cost  per  person  for  fallout  shelter.^-  The  sum  of  the  costs  for 
blast  shelters  and  fallout  shelters  equals  the  total  cost  of  the 
posture  required  to  limit  fatalities  to  an  for  the  given  threat. 
Since  blast  shelters  provide  adequate  fallout  protection,  the  total 
population  has  been  protected  against  fallout.  Therefore,  the  j 
fatalities  bounded  by  a  are  attributed  to  blast  effects. 

E.  ALTERNATIVE  ALLOCATION  PROCEDURE 

Step  4  in  the  present  version  of  the  model  requires  that 
shelters  be  allocated  one  strength  at  a  time,  starting  with  the 
strongest.  This  procedure  shortens  the  allocation  process  because 
the  stronger  shelters  apply  as  assets  when  the  weaker  shelters 
are  being  allocated.  However,  all  shelter  strengths  could  be  con¬ 
sidered  simultaneously  within  the  regions  where  they  are  both 
effective  and  necessary  with  respect  to  each  potential  ground  zero. 
Once  this  was  done  for  the  first  ground  zero,  the  problem  for  suc¬ 
ceeding  ground  zeros  would  be  to  determine  when  and  where  existing 
weaker  shelters  should  be  upgraded  as  a  lower  cost  alternative  .to 
allocating  more  shelters  of  higher  strength  to  previously  unshel¬ 
tered  grid  elements.  The  upgrading  problem  is  avoided  in  the 
present  shelter  allocation  procedure. 

A  situation  is  conceivable  in  which  this  alternative  would  be 
preferred.  In  assigning  relative  priorities  to  determine  which 
grid  elements  will  be  sheltered,  it  is  possible  that  the  cost  per 
person  of  allocating  stronger  shelter  to  one  grid  element  would  be 
found  to  be  lower — due  to  the  applicability  of  more  economical 

1.  See  Appendix  C  for  cost  data  applied  in  pilot  studies. 
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shelter  capacities — than  the  cost  of  allocating  shelter  to  a  less 
densely  populated  element,  where  a  lower  strength  would  provide 
the  required  protection.  This  problem  should  be  considered  in 
further  development  of  the  shelter  allocation  model. 


15 


Ill 


DAMAGE -LIMITING  SHELTER  POSTURES  FOR  HOUSTON: 
A  PILOT  STUDY 


Twenty-five  alternative  shelter  postures  were  developed  and  costed 
for  Houston,  Texas,  to  demonstrate  the  cost/effectiveness  of  a 
damage-limiting  strategy  implemented  through  the  shelter- 
allocation  model  described  in  Section  II.  The  specifications  for 
these  postures  are  listed  below. 

(1)  Houston  1960  at-home  population  (1,225,898) 

A  single  10-Mt  surface  burst 


(2) 

(3) 

(4) 


Blast-shelter  protection  is  assumed  required  inside 
the  7.3-psi  radius  with  respect  to  a  given  ground  zero 

Light  steel-frame  fallout  shelters  are  assumed  to  pro¬ 
vide  blast  protection  against  overpressures  not  greater 
than  7.3  psi.1 2 3 

2 

Alternative  shelter  configurations  and  fatality  levels  con¬ 
sidered  in  this  exercise  are  given  in  Table  1. 


A.  SHELTER  ALLOCATIONS 

Figure  2  shows  the  1960  at-home  population  of  Houston,  Texas, 
allocated  to  the  65  x  65  matrix  of  one-kilometer-square  elements 
used  in  these  studies.^ 


1.  The  overpressure  at  which  the  kill  probability  is  0.5 
for  light  steel- frame  buildings,  according  to  Dikewood’s  blast 
mortality  curves.  L.  Wayne  Davis  et  al.,  Prediction  of  Urban 
Casualties  from  the  Immediate  Effects  of  a~~Nuclear  Attack  (u) 
CONFIDENTIAL,  Dikewood  Corporation,  Contract  No.  OCD-OS-62-203 
(Albuquerque:  April  1963). 

2.  A  listing  of  the  various  shelters  considered  available  for 
these  postures  and  their  respective  costs  is  provided  in  Appendix  C. 

3.  Note  that  the  entire  65  x  65  matrix  is  not  shown  in  any 

of  these  Figures.  Some  unpopulated  rows  at  the  top  and  bottom  have 
been  omitted. 
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Table  1 


SHELTER  POSTURES  CONSIDERED  FOR  HOUSTON 
(Single  10-Mt  Surface  Burst,  At-Home  Population) 


Maximum  Fatality  Level  (a)  Shelter  Posture  Evaluated  to  Achieve  _  ga>^ 


Percent  of 

Number  of 

300 

100 

35 

300/100/ 

300/100 

100/35 

Fallout 

Population 

Fatalities 

psi 

psi 

Psi  , 

35  psi 

psi 

psi 

(7.3  psi 

2.2° 

26,955 

X 

x(3) 

X 

4.0 

49,036 

X 

X 

X 

5.04d 

61,774 

X 

X 

x(4) 

10.0 

122,590 

X 

X 

X 

10.68e 

130,965 

X 

X 

x(5) 

20.0 

245,180 

X 

X 

x(6) 

30.0 

368,769 

X 

X 

X 

35.0 

429,065 

X 

X 

X 

36.67 

449,537 

X 

a.  PF  >  40  fallout  shelters  with  a  blast  rating  of  7.3  psi  are 
allocated  to  all  populations  that  do  not  receive  blast  shelters. 

b.  Numbers  in  parentheses  correspond  to  Figures  depicting  these 
postures . 

c.  the  minimum,  for  300-psi  shelters. 

d.  aj,  the  minimum,  for  100-psi  shelters. 

e.  a^,  the  minimum,  for  35-psi  shelters. 
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Shelter  postures  showing  the  strength  of  the  shelter  allocated 
to  each  grid  element  are  shown  in  Figure  3  through  6  for  four  of 
the  more  interesting  postures  listed  in  Table  1.  The  population 
detail  of  Figure  2  is  not  shown  in  these  Figures,  but  each  matrix 
element  can  be  identified  by  its  column  and  row  coordinates  at  the 
top  and  left  of  each  Figure.  The  first  three  of  these  postures 
provide  the  highest  survival  levels  attainable  with  three  of  the 
different  shelter  configurations  studied:  300/100/35  psi  (Figure  3), 
100/35  psi  (Figure  4)  and  35  psi  (Figure  5).  Figure  6  illustrates 
a  35-psi  posture  designed  (under  assumed  budgetary  constraints)  for 
only  80  percent  survival.  This  posture,  designed  for  a  survival 
level  lower  than  the  highest  attainable,  demonstrates  the  capability 
of  the  shelter  allocation  model  to  accommodate  budgetary  constraints. 
Compare  this  posture  with  the  posture  shown  in  Figure  5,  which  was 
designed  to  achieve  the  highest  survival  level  attainable  with 
35-psi  shelters,  i.e.,  89  percent. 

B.  COST/EFFECTIVENESS 

The  pilot  studies,  based  on  cost  data  described  in  Appendix  B, 
revealed  that  the  35-  and  100-psi  postures  under  the  same  design  con¬ 
ditions  have  substantially  equal  cost/effectiveness  ratios  at  the  65 
to  89  percent  survival  levels.  In  comparison  to  a  35-psi  posture, 
the  100-psi  posture  required  fewer  spaces  for  the  same  level  of 
effectiveness  because  it  offers  a  higher  survival  rate  per  shelter 
space.  This  increased  effectiveness  offsets  the  higher  cost  of  the 
100-psi  shelters.  This  relationship  also  held  for  the  mixed- 
shelter  cases,  i.e.,  the  100-psi  and  300/100-psi  postures  proved  to 
be  as  efficient  as  the  100/35-psi  and  the  300/100/35-psi  postures 
respectively.  This  is  a  significant  finding  because,  as  will  be 
shown  later  (Section  III  C)  100-psi  postures  designed  for  a  single 
10-Mt  attack  prove  far  more  efficient  against  higher-yield  weapons 
than  comparable  35-psi  postures. 
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FIGURE  5  Shelter  Posture:  35-psi  Blast  She 
Maximum  Allowable  Fatalities: 
percent  (One  10-Mt  Surface  Bun 
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FIGURE  6  Shelter  Posture:  35-psi  Blast  Sh« 
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The  35-psi  configuration,  which  shows  no  special  utility  at  the 
65  to  89  percent  survival  levels,  drops  out  of  competition  above  the 
89  percent  level,  the  highest  attainable  when  restricted  to  35-psi 
blast  protection.  At  the  90  to  95  percent  levels,  the  100-psi  pos¬ 
tures  are  more  efficient  than  the  300-  or  300/100 -psi  alternatives. 
The  highest  survival  level  attainable  with  100-psi  shelters  is  95 
percent  under  the  stated  attack  conditions. 

Mixed  postures  comprising  some  300 -psi  shelters  are  the  most 
efficient  at  survival  levels  attainable  above  95  percent.  For 
example,  at  the  98  percent  level  —  the  highest  that  can  be  guaran¬ 
teed  with  300-psi  shelters  --  the  required  300/100 -psi  mix  costs 
only  $548  million  as  compared  to  $828  million  for  a  posture  in  which 
the  choice  of  blast  shelter  is  limited  to  the  300-psi  design.  The 
damage -limiting  mix  compares  even  more  favorably  with  the  universal 
300-psi  posture,  which  costs  twice  as  much  ($1,014  million)  but  which 
could  not  guarantee  more  than  98  percent  minimum  survival  under  the 
same  attack  conditions. 

Figure  7  depicts  the  cost  and  effectiveness  of  the  most 
efficient  postures  developed  in  the  pilot  study.1 2  Postures  are 
included  for  a  number  of  different  survival  levels  ranging  from  65 
to  98  percent.  The  cost/effectiveness  relationship  for  each  of  the 
following  universal  postures  (uniform  protection  for  all)  is  also 
shown  for  comparison: 

(1)  7.3-psi,  PF  =  40  fallout  shelters  (63%  minimum  survival) 

(2)  35-psi  blast  shelters  (89%  minimum  survival) 

(3)  100-psi  blast  shelters  (95%  minimum  survival) 

(4)  300-psi  blast  shelters  (98%  minimum  survival) 


1.  The  numbers  and  costs  of  blast  and  fallout  shelter  spaces 
comprising  each  of  these  postures  are  provided  in  Appendix  D. 

2.  It  is  important  to  recall  that,  as  described  in  Section  II, 
people  who  do  not  receive  blast  shelter  protection  under  the  damage - 
limiting  allocation  procedure  are  allocated  fallout  shelter  protection 
(7.3  psi,  PF  =  40)  at  substantially  lower  cost  per  space  required. 
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FIGURE  7  Cost  and  Effectiveness  of  Alternative  Shelter  Postures 

(single  10-MT  surface  burst  against  the  Houston  at-home  population) 
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As  Indicated  above,  blast  shelters  are  not  required  for  sur¬ 
vival’  levels  less  than  or  equal  to  63  percent.  Fallout  shelters 
alone,  allocated  in  the  grid  elements  Where  the  populations  reside 
and  fully  occupied,^  would  enable  63  percent  survival.  Although  not 
shown  in  Figure  7,  only  18  percent  of  the  Houston  population  would 
survive  if  no  special  shelter  were  provided  and  people  remained  in 
their  residences. 

Average  and  marginal  costs  per  survivor  for  each  shelter  con¬ 
figuration  are  readily  obtainable  from  the  model  output. 

C.  YIELD  SENSITIVITY 

Two  alternative  attacks  were  used  to  test  the  yield  sensitivity 
of  several  of  the  postures  developed  for  this  pilot  study.  Designed 
to  counter  a  single  10 -Mt  threat  to  the  Houston  at-home  population, 
the  postures  were  tested  against  single  3-Mt  and  100 -Mt  weapons, 
each  a  surface  burst  at  its  optimum  ground  zero.  The  results  for 
the  best  10 -Mt  postures  —  those  which  were  determined  to  be  most 
efficient  at  each  survival  level  studied  (Figure  7)  —  are  shown  in 
Figure  8. 

As  reported  earlier,  the  efficiency  of  alternative  postures 
comprising  35-psi  and  100 -psi  blast  shelters  was  determined  to  be 
substantially  equivalent  against  the  design  threat  of  10  Mt.  However, 
the  efficiency  of  35-psi  postures  is  significantly  below  that  of  the 
100-psi  postures  when  subjected  to  the  higher  yield  (100-Mt)  attack. 
The  yield  sensitivity  of  35-  and  100-psi  shelter  designs,  having 
comparable  10 -Mt  survival  levels  at  each  of  four  budget  levels,  is 
reflected  in  Table  2. 


1.  Occupancy  percentages  can  be  varied  in  applying  the  model. 
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Table  2 


YIELD  SENSITIVITY  OF  COMPARABIE  35-PSI  AND 
100-PSI  BIAST  POSTURES 

(Designed  to  Counter  a  Single  10-Mt  Surface  Burst, 
Houston  At-Home  Population) 


Budget  Level 
millions  of  $ 


Blast  Shelter 
Utilized  |n 
Posture" 


Minimum  Survival  Level,  % 
On-Design  Of^-Design 

(10-Mt)  3-Mt  100 -Mt 


290 

35-psi 

89 

93 

38 

100-psi 

89 

93 

47 

200 

35-psi 

80 

88 

22 

100-psi 

80 

89 

29 

145 

35-psi 

70 

83 

18 

100-psi 

70 

83 

23 

120 

35-psi 

65 

82 

16 

100-psi 

65 

82 

19 
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People  for  whom  blast  shelter  is  not  provided  are 
allocated  fallout  shelter  protection:  7.3  psi,  PF=40 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  pilot  study  reported  in  Section  III  confirms  that  shelter 
postures  developed  under  the  damage-limiting  strategy  proposed 
here  can  make  it  possible  to  limit  fatalities  to  a  stipulated  level 
regardless  of  where  the  weapon  may  be  delivered  within  a  target  city. 
This  study  also  established  that  by  allocating  shelters  according 
to  the  model  presented,  fatalities  can  be  limited  to  the  level 
expected  if  shelters  of  the  highest  psi  rating  under  consideration 
were  provided  universally  (at  much  greater  cost)  and  if  the  weapon 
were  detonated  at  the  enemy’s  best  aimpoint. 

It  is  recommended  that  the  model  be  used  to  test  the  sensitivity 
of  the  results  to  all  input  criteria  and  parameters  which  signifi¬ 
cantly  affect  shelter  requirements  and  costs.  Initial  plans  are  that 
the  model  will  be  used  to  study  the  sensitivity  of  the  results  to 
the  size  of  the  matrix  elements  (now  one-kilometer  square).  Alter¬ 
natives  include  the  use  of  two-kilometer-square  elements,  which 
would  reduce  the  total  number  of  matrix  cells  from  4,225  to  1,024; 
or  the  use  of  census  tracts,  which  are  not  uniform  in  size  or  shape. 
The  geographic  center  of  each  tract  can  be  used  as  the  basis  for 
measuring  weapon  effects  over  the  whole  tract  area. 

The  model  also  will  be  used  to  develop  protective  postures  for 
different  population  modes:  at-home,  at-work  or  rush-hour.  These 
postures  will  then  be  tested  against  alternative  times  of  attack : 
day,  night  or  rush  hour,  as  well  as  alternative  yields  and  numbers 
of  weapons.  The  results  will  delineate  those  postures  with  the 
lowest  cost/effectiveness  ratio  over  a  range  of  circumscances . 

Some  possible  avenues  of  research  for  further  development  of 

the  model  itself1  are  as  follows: 

—  . - .  — .  *■ 

1.  Also  see  discussion  in  Section  HE. 
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(1)  Incorporate  a  probability  of  kill  (P,  )  function  for  use  in 
estimating  fatalities  among  people  occupying  different  strength 
shelters.  The  P.  would  vary  according  to  the  specific  over¬ 
pressure  level  applied  and  should,  of  course,  be  consistent  with 
the  design  specifications  for  the  shelter.  A  cookie-cut  criteri< 
is  applied  in  the  present  model,  i.e.,  all  people  live  if  the 
overpressure  does  not  exceed  the  psi  rating  of  their  shelter;  al] 
people  die  if  it  does. 

(2)  Allow  for  divisibility  or  combination  of  grid  elements  in 
the  assignment  of  shelter.  In  the  present  model,  the  total 
population  of  a  grid  element  is  the  basic  decision  unit  and 
receives  a  uniform  level  of  protection.  This  procedure  does  not 
generally  permit  the  allocation  of  shelter  in  the  most  economica] 
sequence  or  priority,  given  the  existence  of  budgetary  constraint 
For  example,  if  grid  elements  were  not  treated  as  indivisible, 
shelters  of  the  most  economical  size,  e.g.,  1,000-man  occupancy, 
could  be  utilized  to  the  maximum  extent  possible  before  allocatir 
any  of  the  smaller  shelters  which  are  now  necessary  to  provide 
full  protection  for  all  the  people  in  an  element.  Alternatively, 
sparsely  populated  elements  could  be  joined  together  for  protec¬ 
tion  purposes  in  order  to  permit  more  efficient  shelter  assignmer 

Clearly,  a  program  that  would  allocate  different  protection  level 
to  different  populations  within  the  same  area  has  important  social  and 
political  implications  which  must  be  understood  and  accommodated  befor 
shelters  could  be  constructed  and  installed.  Equitable  and  politicall 
acceptable  allocation  is  only  one  problem.  Feasible  shelter  sites 
would  have  to  be  identified  through  engineering  surveys;  and  property 
ownership  problems,  including  rights-of-way,  would  have  to  be  resolved 
before  sites  could  be  selected.  A  model  such  as  that  presented  in  thi 
Study  would  be  an  invaluable  tool  for  evaluating  alternative  courses 
and  decisions  at  each  point  in  the  development  of  this  type  of  operat¬ 
ing  program. 
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*  M 


GLOSSARY  OF  NOTATION 


a..:  Matrix  element  in  the  i-th  row,  j-th  column 

being  evaluated  for  vulnerability  or  sheltering 

a  :  Matrix  element  in  the  m-th  row,  n-th  column 

being  evaluated  as  an  Actual  Ground  Zero  (AGZ) 


ci-ik:  Avera?e  cost  Per  ‘''Pace  for  protecting  the  popu 
-1  lation  in  the  i-,  j-th  element  with  k  strength 
shelter.  This  average  cost  reflects  the  most 
economical  combination  of  shelter  sizes  that 
will  accommodate  that  number  of  people. 

d .  . :  Distance  of  an  element,  a. .,  from  the  target 

J  element,  a  1-5 

mn 

d..  =  V(i-m)5  +  (j-n)2 


The  allowable  fatalities  in  R^ 


P 


o 


£  p . . ;  The  sum  of  the  population  in  the  unprotectable 
RQ  J  region,  RQ 


Pij  :  The  sum  of  the  population  within  r^ 

> 

p.  Population  in  element  a.. 

13  13 

R^:  A  region  within  which  k  shelter  strength  will 
A  protect  and  a  lower  psi  shelter  will  not 

(k  =  0,  1,  2,  ...»  N)  in  relation  to  a  specific 
AGZ 

Rq:  Region  about  AGZ  within  which  strongest  shelter 
under  consideration  (S2)  will  not  provide  ade¬ 
quate  blast  protection  against  a  given  yield 
weapon.  Rq  contains  all  matrix  elements  i,j 
located  at  distances  from  ground  zero  where 
the  blast  overpressure  is  greater  than  the  psi 
rating  of  the  strongest  shelter  under  consideration 
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Outer  limit  of  R.  as  a  distance  from  ground 
zero 


r  :  Outer  limit  of  R  as  a  distance  from  ground 
zero  * 


s,  :  Shelter  strength  being  allocated  (k  =  1  for  the 
strongest  shelter,  2  for  the  next  strongest,  etc 

tlv  =  Z  u  •  •  i  The  total  unprotected  population  within  region  R 
*k  13 

u .  . :  Population  p-jj  of  an  element  a^j  which  does  not 
3  have  adequate  blast  protect  ion,  e.g. ,  the  popu¬ 
lation  of  an  element  in  R2  which  does  not  have 
Sj  or  stronger  (s^)  shelter 

V  :  The  vulnerable  population  with  respect  to 

a  given  ground  zero,  a^,  that  is  unprotectable 
with  the  shelter  strength,  s^,  being  allocated 

Vmn  =  Po  +  U1  +  U2  +  *  •  *  +  \-l 

an:  The  level  at  which  fatalities  are  to  be  upper 
bounded  (n  =  1,  2,  ...,  m) 

a^:  The  minimum  level  at  which  fatalities 
can  be  upper- bounded  with  the  strongest 
shelter,  s^,  under  consideration 

a^:  The  maximum  fatality  level  for  which  at 
least  some  of  the  strongest  shelters,  s,, 
would  be  required.  1 
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Appendix  B 


SPECIFICATIONS  FOR  PROGRAM  DAMIT1 

This  appendix  describes  the  shelter  allocation  model  in  sufficient 
detail  to  provide  a  more  thorough  understanding  of  the  basic  process 
and  to  facilitate  computer  programming.  Although  these  specifica¬ 
tions  incorporate  refinements  to  the  initial  computer  program 
(ALLOCATE)  used  for  the  studies  presented  in  Section  III,  the  basic 
allocation  procedures  are  identical. 

The  main  flow  of  Program  DAMIT  is  presented  in  Figure  B-l. 

The  remainder  of  this  appendix  is  devoted  to  a  detailed  description 
of  the  steps  depicted  in  this  diagram. 

I.  INPUT  ROUTINE  (Figure  B-2) 

I-A.  Population  Matrix 

A  magnetic  tape  containing  the  population  matrix  should  be 
prepared.  All  pertinent  dimension  statements  should  be  modified 
to  conform  to  the  size  of  the  population  mctrix.  Total  population 
should  be  computed  immediately  and  stored  for  later  use. 

I-B.  Data  Card(s)  1 

These  are  cards  containing  alphamerical  descriptions  of  the 
problem.  Example:  "1960  HOUSTON  NIGHTTIME  POPULATION."  This  can 
be  simplified  by  assigning  problem  identification  numbers. 

I-C.  Data  Card  2 

Columns  1-5  :  Number  of  rows  in  population  matrix. 

6-10:  Number  of  columns  in  population  matrix. 

11-15:  Number  of  shelter  strengths  in  posture. 


1.  DAMage  limITing  strategy. 
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'  Read  In 
Population 
v  Matrix 


_ J 

Sum 

Population 


Read 
Data  Card(s)  #1 
^[Heading  Card(s)] 


C.  f 
Read 


/  Read  \ 
\Pata  Card  tV 


D. 


Read 

Data  Card(s)  #3 
[Shelter  Information] 


Read 

(Print  Option 
Card 


F. 


FIGURE  B-2  Input  Routine 
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16-20:  Weapon  yield  (megatons). 

21-25:  Prescribed  fatality  level  (a^)  in  percent. 

26-30:  Blast  protection  level  (psi)  given  to  fallout 
shelters . 

31-40:  Cost  per  person  of  fallout  shelters. 

I-D.  Data  Card  3 — Shelter  Information  Card(s) 

For  maximum  input  flexibility,  there  should  be  one  card  for 
each  shelter  size. 

Columns  1-5  :  Shelter  strength  (psi). 

6-10:  Shelter  size  (number  of  persons  capacity). 
11-20:  Shelter  cost  (nearest  dollar). 

I-E.  Print  Option  Card 

This  program  should  permit  the  examination  of  output  at 
critical  stages  of  the  various  iterative  processes.  This  will 
maximize  the  research  potential  of  the  program. 

Each  column  of  the  print  option  card  represents  a  particular 
table  in  the  program.  A  "0"  in  a  column  indicates  that  the  table 
is  not  to  be  printed.  A  "1"  in  the  column  specifies  that  it  is 
to  be  printed. 

Certain  tables  such  as  the  summary  of  inputs  and  the  final 
shelter  posture  should  be  printed  automatically. 

I-F.  Table  1.00 — Summary  of  Inputs 

This  table  indicates  all  information  on  data  cards  1,  2  and  3. 

I-G.  Table  1.01 — Population  Matrix  (optional) 

This  table  is  a  printout  of  the  population  matrix.  The  total 
population  should  be  printed  beneath  the  table  heading.'*' 


1.  See  Subroutine  PRINT1  in  the  Civil  Defense  Project  Subroutine 
Library  Manual.  This  publication  will  henceforth  be  abbreviated  as 
CDPSI/I. 
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II.  DETERMINATION  OF  PSI  RADII  (Figure  B-3)2 

This  portion  of  the  program  determines  the  distance  from 
a  given  ground  zero  that  a  shelter  of  a  certain  strength  will  be 
effective.  These  are  the  r^,  i.e.,  the  outer  limits  of  the  region 
within  which  a  given  strength  shelter  in  the  posture  will  protect 
and  a  lower  strength  shelter  in  the  posture  will  not. 

II-A.  Prepare  Psi  Array 

The  shelter  strengths  read  in  at  Steps  I-C  and  I-D  are  collected 
into  an  array  in  descending  order  of  magnitude.  The  number  of  radii 
to  be  computed  must  also  be  determined. 

II-B.  Call  Subroutine  RADIUS 

A  description  of  this  subroutine  may  be  found  in  the  CDPSU4. 

The  arguments  for  this  subroutine  are 

(1)  Number  of  radii  required 

(2)  Array  of  radii  (output) 

(3)  Psi  array  (input) 

(4)  Yield  (megatons) 

II-C.  Store  Radii 

The  array  of  radii  should  be  stored  for  later  used  in  the 
program. 

II-D.  Print  (optionally)  psi  and  radii. 

III.  TEST  FOR  CONSISTENCY  OF  INPUT  (Figure  B-l)2 

The  fatality  level,  an>  inputted  via  data  card  2  must  be  con¬ 
sistent  with  the  shelter  configuration  prescribed  on  the  shelter 
information  cards.  If  is  so  large  that  the  highest  shelter 
strength  inputted  would  not  be  allocated,  a  great  deal  of  computer 


1.  This  represents  a  modification  of  the  operational  version 
of  the  Program  (ALLOCATE)  where  the  radii  are  required  as  input. 

2.  This  represents  a  modification  of  Program  ALLOCATE, 
where  no  such  test  is  made. 
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FIGURE  B-3  Determination  of 


time  would  be  wasted  in  trying  to  allocate  the  highest  strength 
shelters  when  such  shelters  are  actually  not  needed.  This  test 
eliminates  this  possibility  by  printing  out  an  error  message  and 
halting  the  operation  if  the  specified  fatality  level  is  equal  to 
or  greater  than  the  upper  bound  of  fatalities  associated  with  the 
next  highest  strength  shelter. 

III-A.  Determine  Matrix 

is  the  population  in  the  region  around  a  given  ground  zero 
within  which  the  next-to-strongest  shelter  in  the  posture  will  not 
provide  adequate  blast  protection  against  a  given  yield  weapon. 

The  P^  matrix  represents  the  set  of  P^  values  determined  when  each 
element  in  the  population  matrix  is  assumed  to  be  the  target. 

For  example,  Table  B-l  presents  a  hypothetical  5x5  population 
matrix.  One  kilometer  separates  each  adjacent  element.  Suppose 
that  element  (3,3)  is  regarded  as  the  ground  zero  and  that  the 
distance  within  which  no  protection  can  be  given  by  the  next-to- 
highest  strength  shelter  available  is  one  kilometer.  This  means 
that  the  population  in  elements  (2,3),  (3,2),  (3,3),  (3,4),  and 
(4,3)  cannot  be  protected  P1(3,3)  =  21.  All  other  elements  are 
at  a  distance  greater  than  one  kilometer  from  ground  zero  (3,3). 

Pj  (3,3)  is  entered  into  Table  B-2  and  circled. 


Table  B-l 

Population  Matrix 

X 

1  2  3  4  5 

1 

1112  1 

2 

2  13  3  1 

3 

0  5  ©  4  3 

4 

1  4  4  4  2 

5 

2  3  2  2  1 

Now  suppose  that  element  (1,1)  is  the  ground  zero.  The  popu¬ 
lation  one  kilometer  or  less  from  this  target  is  represented  by 
elements  (1,1),  (1,2),  and  (2,1).  Therefore,  P1  (1,1)  =  4,  and  this 
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value  is  circled  in  Table  B-2.  In  the  same  fashion,  all  other  ele¬ 
ments  are  assumed  to  be  ground  zeros  and  the  values  of  (I,J) 
determined.  Subroutine  PSUBK  is  designed  to  obtain  the  P^  matrix. 
More  generally,  it  is  designed  to  obtain  any  such  matrix,  P^,  where 
the  population  within  a  given  radius  from  each  possible  ground  zero 
is  to  be  determined.1 2 


Table  B-2 
P1  Matrix 


X 

1 

2 

3 

4 

5 

1 

© 

4 

7 

7 

4 

2 

4 

12 

13 

13 

8 

3 

8 

15 

@ 

19 

10 

4 

7 

17 

19 

16 

10 

5 

6 

11 

11 

9 

5 

III-B.  Determine  P^CMAX) 

Subroutine  PKMAX  uses  the  P^  matrix  as  input  and  determines 
the  maximum  value  in  the  matrix. ^  In  the  above  example,  P^CMAX)  =  21, 
and  this  would  be  the  output  of  the  routine. 

III-C.  Test  if  an  Is  Greater  than  or  Equal  to  P^(MAX) 

If  the  inputted  fatality  level  aR  is  greater  than  or  equal  to 
P^(MAX),  then  the  program  would  unnecessarily  attempt  to  allocate 
the  highest  strength  shelter  and  computer  time  would  be  wasted. 

To  prevent  this,  the  following  error  message  should  be  printed 
(off  line)  if  this  condition  exists: 

”  Pl(MAX)  =  _ .  INPUTTED  FATALITY  LEVEL  =  _  WHICH  IS  GREATER 

THAN  OR  EQUAL  TO  Pl(MAX).  NO  SI  SHELTERS  CAN  BE  ALLOCATED. 

CHANGE  EITHER  ALPHA  N  LEVEL  OR  STRENGTH  OF  SI  SHELTER. " 

The  program  should  be  terminated  after  this  message  is  printed. 

1.  See  Subroutine  PSUBK  in  the  CDPSLM . 

2.  See  Subroutine  PKMAX  in  the  CDPSIi-1. 
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IV.  SET  UP  INITIAL  CONDITIONS  FOR  LOOP  ENTRY 

The  allocation  begins  at  the  element  with  the  maximum  unprotect- 
able  population  level.  In  this  case,  the  VQ  matrix  is  identical  to 
the  PQ  matrix  and  VQ(MAX)  =  Pq(MAX)  =  ar 

IV-A.  Determine  V  Matrix 
-  o  - 

Here  Subroutine  PSUBK  can  be  called  again  using  as  input  the 
radius  within  which  no  person  can  be  protected  by  the  highest 
strength  shelter.'*' 

IV-B .  Determine 

Subroutine  PKMAX  may  be  used  to  obtain  VQ(MAX)  =  a^.  The 
coordinate  of  this  maximum  must  be  stored  for  use  in  the  allocation 
loop,  should  be  stored  for  final  output. 

IV-C.  Test  if  a  Is  Less  than  a, 

— . .  n  — . .  ■  - . — —  1 

If  the  inputted  fatality  level,  an,  is  less  than  a^,  then  it 
is  assumed  that  is  the  fatality  level  prescribed,  for  it  is 
impossible  to  guarantee  fatalities  below  the  level. 

IV- D.  Initialize  Vq(MAX)  and  K 

Vq(MAX)  is  set  to  an  and  K  is  set  equal  to  1. 

V.  ALLOCATE  SHELTERS  OF  STRENGTH  sk 

Part  V  represents  the  shelter  allocation  loop.  The  individual 
steps  are  described  in  detail  below. 

V- A.  Determine  fh(M,M)  and  U?(M,N) 

Subroutine  FSUBK  (see  Figure  B-4)  is  designed  to  obtain  for 
a  specific  target  (1)  fk  =  the  number  of  persons  that  may  be  left 
unsheltered  without  exceeding  the  fatality  limit,  (2)  Uk  =  the  number 
of  persons  in  the  protectable  region  that  have  not  been  given 
adequate  shelters,  and  (3)  Vk  =  PQ  +  U^,  where  PQ  is  the  population 
that  cannot  be  protected  even  with  the  highest  strength  shelter. 

Vk  may  be  thought  of  as  the  "vulnerable”  population. 

1.  The  P1  matrix  is  no  longer  needed.  Store  VQ  in  the  same  area. 
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The  arguments  to  the  subroutine  are  as  follows : 

Inputs 

M:  Row  coordinate  of  the  ground  zero  of  interest 

N:  Column  coordinate  of  the  ground  zero  of  interest 
NR:  Number  of  rows  in  population  matrix 
NC:  Number  of  columns  in  population  matrix 

Po(M,N) :  Unprotectable  population  for  ground  zero  of  interest 
NL:  Number  of  radii  to  be  considered  (=  K  +  1) 

RAD:  Array  of  radii  computed  in  Step  II 
SHELT:  Array  of  assigned  shelter  strengths 
POP:  Population  array 
SK:  Strength  of  shelter  being  allocated 
AN:  Maximum  allowable  fatalities,  a 

n 

Outputs 

F:  f^  -  allowable  fatalities 

U:  =  unprotected  population  in  protectable  region 

V:  =  vulnerable  population 

The  basic  steps  of  the  routine  are: 

Step  1:  Reduce  Matrix  to  Minimum  Needed.  Nearly  all  elements 
whose  distances  are  greater  than  the  largest  inputted  radius, 
RAD(NJ),  are  eliminated  from  consideration.  Precautions  are  taken 
to  ensure  that  the  limits  of  the  reduced  matrix  are  within  the 
population  matrix. 

Step  2:  Clear  SUM  (accumulator  for  computation). 

Step  3:  L  =  1.  Thf  first  radius  to  consider  is  RAD(l)  =  rQ. 

L  is  the  radius  index.  We  are  now  in  the  unprotected  population 
accumulation  loop,  starting  at  the  upper  left-hand  coordinates  of 
the  reduced  matrix. 

Step  4:  Determine  the  Distance  of  the  Center  of  the  Element 

from  the  Ground  Zero.  _ 

D(I, J)  =  I  -  M)*  +  (J  -  N)2 
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St *'p  *ii  Test  if  D(I,J)  >  RAD(L)  =  RAD(L)  =  rQ.  If  this  dis¬ 
tance  is  **4ual  to  or  less  than  RAD(l),  then  the  element  is  in  the 
unprotectable  region,  and  its  population  would  be  included  in  PQ(M,N). 
If  it  is  greater  than  RAD(l),  the  next  test  must  be  made. 

Step  6:  Test  if  D(I,J)  >  RAD(L  +  1).  At  the  first  iteration, 
the  test  determines  whether  D(I,J)  >  RAD(2)  =  r^.  If  it  is,  then  it 
tests  whether  D(I,J)  >  RAD(3)  =  r2>  and  so  on  until  the  largest  radius, 
RAD(NL)  to  be  considered  has  been  reached.  If  the  distance  is  greater 
than  RAD(NL) ,  the  element  is  outside  the  region  of  interest,  and  the 
program  branches  to 
interest,  control  go 

Step  7:  Test  if  Element  is  Adequately  Sneltered.  Using  the 
array  of  assigned  shelters,  test  to  determine  whether  a  shelter 
strength  that  may  have  been  assigned  to  the  element  when  other 
targets  were  considered  is  strong  enough  to  withstand  the  psi  level 
it  will  encounter  for  this  ground  zero.  If  it  is,  the  population 
in  the  element  is  adequately  sheltered  for  the  ground  zero  under 
consideration.  If  it  is  not,  control  goes  to  Step  8. 


Eep  9.  If  the  element  is  within  the  region  of 
,  to  Step  7. 


Step  8:  Accumulate  Unprotected  Population. 


Step  9.  Test  if  More  Elements  Are  to  be  Considered.  If  yes, 
L  is  restored  to  equajL  1,  and  a  new  element  in  the  reduced  matrix 
is  considered.  If  no,  go  to  Step  10. 


V-B. 


Step  10:  U^^(M,N)  =  the  Sum  of  the  Unprotected  Population,  SUM. 

Step  11:  Vj;(M,N)  =  the  Vulnerable  Population  =  PQ(M,N)  +  U^(M,N), 
Step  12:  f^(M,N)  =  Allowable  Fatalities  =  g  -  V^(M,N). 

Test  if  Uk(M,N)  >  fk(M,N) 

If  the  unprotected  population  in  the  protectable  region  is 


greater  than  the  allowable  fatalities,  then  enough  s^  shelters  must 

be  assigned  to  keep  fatalities  within  the  bounds  of  an  and  the 

then  no  s,  shelters 
k 

need  to  be  allocated  to  keep  fatalities  within  the  bounds  of  a  ; 

n 

therefore,  the  program  skips  to  step  V-D. 


program  goes  to  step  V-C.  If  U^(M,N)  <  f^(M,N) 
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V-C.  Allocate  Shelters 

Subroutine  SHELTER  (see  Figure  B-5)  is  designed  to  allocate 
shelters  of  a  given  strength  and  given  capacities  to  the  population 
of  a  region  around  a  prescribed  ground  zero  where  the  shelter 
strength  is  adequate  to  protect  the  population.  The  routine  assigns 
shelters  only  up  to  a  specified  population  level  within  the  region 
and  assignments  are  made  to  those  elements  where  the  average  cost 
per  pens  is  smallest. 1  Subroutine  MINLCOST  is  used  to  determine 
the  average  cost  per  person  of  each  element  in  the  region. 

The  arguments  far  Subroutine  SHELTER  are  as  follows ; 

!  '  Inputs 


II 

JJ 

R1 

R2 

NC 

NR 

IALF 

JSS 

NSM 

NSS 

JS1Z 

JCOST 

IPOP 

IWORK1 

IW0RX2 

IW0RK3 

IWOEK4 

IW0RK5 

IW0RK6 

NRNC 


Row  coordinates  of  ground  zero 

Column  coordinates  of  ground  zero 

Inner  radius  of  region  of  interest  (kilometers) 

Outer  radius  of  region  of  interest  (kilometers) 

Number  of  columns  in  population  matrix 

Number  of  rows  in  population  matrix 

Allowable  fatalities  for  this  region 

Shelter  strength  of  shelters  being  allocated 

Smallest  size  shelter  in  the  shelter  size  array 

Number  of  shelters  in  the  shelter  size  array 

Shelter  size  array 

Shelter  cost  array 

Population  array 

Dummy  work  area 

Dummy  work  area 

Dummy  work  area 

Dummy  work  area 

Dummy  work  area 

Dummy  work  area 

NR  times  NC 


1.  This  represents  a  modification  of  Program  ALLOCATE, 
where  the  average  cost  per  required  space  is  utilized. 
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Output 

MSHELT:  Shelter  strength  array 

Figure  B-5  indicates  the  procedure. 

Step  1:  All  work  arrays  are  cleared  and  variables  initialized 
where  necessary. 

Step  2:  The  two  inputted  radii  are  squared. 

Step  3:  The  size  of  the  population  matrix  is  reduced  to 
approximately  those  elements  whose  distances  are  equal  to  or  within 
the  outer  radius,  r^,  from  the  inputted  ground  zero.  Precautions 
are  taken  to  ensure  that  the  limits  of  the  reduced  matrix  are 
within  the  population  matrix. 

Step  4:  If  the  square  of  the  distance  of  the  center  of  the 

2  2 

element  is  greater  than  r^  and  less  than  or  equal  to  r  ,  the 
element  is  in  the  region  of  interest  and  the  next  instruction  is 
at  Step  5.  Otherwise,  the  next  element  is  considered. 

Step  5:  The  element  is  tested  to  determine  whether  it  has 
already  been  assigned  adequate  shelter.  If  not,  the  next  instruc¬ 
tion  is  at  Step  6.  Otherwise,  the  next  element  is  considered. 

Step  6 :  The  population  and  coordinates  are  stored  in  two 
corresponding  arrays. 

Step  7:  The  minimum  average  cost  per  person  is  determined 
using  Subroutine  MINLCOST.  This  routine  (see  Figure  B-6)  is 
designed  to  determine  the  average  cost  per  person  of  sheltering 
the  population  of  an  element  using  that  combination  of  shelter 
sizes  and  costs  that  will  yield  the  minimum  cost  per  person.  The 
routine  also  returns  the  shelter  size  distribution  that  yields  the 
minimum  cost  per  person. 

The  arguments  for  Subroutine  MINLCOST  are  as  follows : 

Inputs 

NPOP:  Population  to  be  sheltered 
NSM:  Smallest  size  shelter  in  the  array  of  shelter  sizes 
NSS:  Number  of  shelter  sizes  in  the  array  of  shelter  sizes 
NS1Z:  Shelter  size  array 
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FIGURE  B-6  Subroutine  Minlcost 
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NCOST:  Shelter  cost  array 
IWORK4:  Dummy  work  array 
JA:  Dummy  work  array 
IW0RK7 :  Dummy  work  array 


Outputs 


IS1ZE:  Array  giving  number  of  shelters  of  each  size  that 
yields  minimum  average  cost 

MIN:  Minimum  cost  per  person 


Figure  B-6  indicates  the  procedure 

(1)  All  work  and  output  areas  are  cleared. 

(2)  The  initial  size  distribution  is  determined  by: 

(a)  Assigning  the  largest  size  shelter  until  popu¬ 
lation  is  exceeded; 

(b)  Reducing  the  number  of  largest  size  shelters  by 
one  and  assigning  the  next  largest  to  the  remaining 
population  until  this  remainder  is  exceeded; 

(c)  Reducing  the  number  of  the  next-to-largest 
shelters  by  one  and  assigning  the  next  size  shelter 
to  the  remainder,  and  so  on. 

(3)  The  total  cost  of  the  resulting  shelter  size  distri¬ 
bution  is  determined  and  divided  by  the  size  of  population. 

(4)  The  average  cost  of  using  a  larger  shelter  instead  of 
several  smaller  ones  is  determined.  For  example,  one 
500-man  shelter  may  be  used  instead  of  four  100-man  shelters 

(5)  This  is  compared  with  the  initial  average.  If  the 
alternative  cost  is  smaller,  the  initial  shelter  distri¬ 
bution  is  replaced  by  the  new  distribution.  Subsequent 
costs  and  distributions  are  compared  with  the  results  of 
previous  iterations. 

(6)  The  subroutine  returns  the  minimum  cost  determined 
in  (5)  and  the  final  size  distribution  yielding  this  cost. 

Step  8;  This  cost  is  stored  in  an  array. 

Step  9:  If  more  elements  are  i~  the  matrix,  Steps  4  through  8 
are  repeated.  Otherwise,  the  routine  goes  to  Step  10. 

Step  10:  The  cost  array  is  searched  for  the  highest  cost 
element . 

Step  11:  If  two  or  more  elements  have  equally  high  costs,  the 
one  with  the  smallest  population  is  chosen. 
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Step  12:  The  highest  cost  elements  are  accumulated,  until 
Step  13:  The  sum  exceeds  the  allowable  fatalities,  IALF. 

Step  14:  The  population  of  the  last  element  added  is  subtracted 
from  this  sum. 

Step  15:  The  population  in  the  remaining  elements  in  the  popu¬ 
lation  array  is  given  shelter  protection  of  strength  JSS.  This  is 
stored  in  the  array  called  MSHELT,  and  the  routine  returns  to  the 
calling  program. 

V-D.  Test  if  There  Are  More  Ground  Zeros  for  Given  s^ 

If  all  possible  ground  zeros  have  not  been  considered,  the 
ground  zero  with  the  next  highest  vulnerability  in  the  matrix 
(as  determined  at  either  Step  IV-A  or  VI-B)  is  chosen.  This  is 
accomplished  by  first  setting  the  old  maximum  to  zero  and  then 
using  Subroutine  PKMAX  to  determine  the  coordinates  of  the  new 
maximum. 

V-E.  Test  if  There  Are  More  Shelters  to  Allocate 

If  all  of  the  snelter  strengths  have  not  been  allocated,  the 
program  goes  to  Step  VI.  If  all  have  been  allocated,  the  program 
goes  to  Step  VII. 

VI.  DETERMINATION  OF  NEW  MATRIX 

When  the  next  strength  shelter  is  to  be  allocated,  a  new 
vulnerability  matrix  must  be  obtained.  This  is  accomplished  by 
successive  application  of  Subroutine  FSUBK  (Step  VI-B)  and 
replacing  the  old  V^^  values  with  the  new  V^  values  (Step  VI-C). 

When  all  of  the  V^  values  have  been  computed,  the  coordinates  of 
the  maximum  value  of  the  V^  matrix  are  obtained.  This  is  the  first 
ground  zero  for  the  new  shelter  allocation.  Control  returns  to 
Step  V. 
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VII.  GIVE  MINIMUM  BLAST  AND  FALLOUT  PROTECTION  TO  UNPROTECTED 
ELEMENTS 

Minimum  blast  and  fallout  protection  are  given  to  all  elements 
not  receiving  higher  than  minimum  shelters.  The  level  of  the 
minimum  is  specified  in  data  card  2. 

VIII.  FINAL  OUTPUT  ROUTINE  (Figure  B-7) 

VIIi-A.  Print  Final  Shelter  Strength  Posture  (optionally) 

VIII-B.  Write  Final  Shelter  Matrix  on  Tape 

The  final  shelter  posture  is  written  on  tape  with  proper 
labeling. 

VIII-C.  Print  Details  for  Each  Element  (optionally) 

If  the  option  is  exercised,  the  indicated  list  is  printed. 

In  any  case,  the  details  are  written  on  tape. 

VIII-D,  Print  Final  Summary  of  Posture 

This  is  not  an  optional  step.  It  is  assumed  that  this  will 
always  be  desired  information. 

VIII- E.  Write  Summary  Information  on  Tape 
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^inal  Shelter 
Posture ■/* 
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Final  Shelter 
Matrix.  Use 
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FIGURE  B-7  Final  Output  Routine 


Appendix  C 
SHELTER  COSTS 
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SHELTER  COSTS 


The  cost  data  used  in  the  development  and  costing  of  alternative 
shelter  postures  for  Houston,  Texas,  are  as  follows: 

Cost 

Blast  Shelters  Total  Cost  Per  Space 


35  psi 

' 

100-man 

$  49,557 

$  496 

500-man 

165,624 

331 

1000-man 

324,576 

325 

100  psi 

100-man 

58,060 

581 

500-man 

195,055 

390 

1000-man 

379,147 

379 

300  psi 

100-man 

77,067 
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Cl.  FALLOUT  SHELTERS 

The  657,745  additional  spaces  required  under  a  full  fallout 
shelter  program  for  Houston  are  estimated  to  cost  an  average  of 
$05  per  space.  The  OSA(CD)  minimum  criteria  for  fallout  shelters 
is  PF  =  40.  These  spaces  would  be  located  in  buildings,  specifi¬ 
cally  hospitals  and  schools,  which  are  of  light  steel-frame 
construction  in  Houston. 

C2.  DATA  SOURCES 

The  above  shelter  costs  were  developed  by  the  IDA  Civil 
Defense  Project  Group  in  1964  for  use  in  local  area  studies  focused 
on  Houston.  The  detailed  data  base  including  specific  sources  is 
available  but  as  yet  unpublished.  Specific  queries  in  this  regard 


may  be  directed  to  the  Head,  Civil  Defense  Project,  Economic  and 
Political  Studies  Division,  Institute  for  Defense  Analyses, 

400  Army-Navy  Drive,  Arlington,  Virginia. 

C3.  COST  PER  REQUIRED  SPACE 

The  element -by-element  differences  in  cost  per  required  space 
are  illustrated  below  for  four  elements  requiring  100-psi  protection 
in  reference  to  a  particular  ground  zero. 


Population 

Required 

Number 

Spaces3 

Most  Efficient 
Shelter  Mix 

Shelter  Cost 

Relative 
Priority 
for  Sheltei 
Allocation* 

Unit 

Total 

Per 

Space 

No. 

Occupancy 

810 

900 

1 

l,000c 

$379,147 

$379,147 

$421 

(3) 

995 

1,000 

1 

1,000 

379,147 

379,147 

379 

(1) 

1,130 

1,200 

1 

1,000 

379,147 

379,147 

plus 

2 

100 

58,060 

116,120 

$495,267 

413 

(2) 

700 

700 

1 

500 

195,055 

195,055 

plus 

2 

100 

58,060 

116,120 

$311,175 

444 

(4) 

a.  Populations 

are 

rounded  up 

to  the  next 

hundred 

for  this 

purpose 

because  all  shelter  occupancies  are  in  hundreds. 

b.  Should  a  choice  exist  between  two  elements  having  the  same  cost 
per  required  shelter  space,  the  model  gives  priority  for  shel¬ 
tering  to  the  one  of  highest  population  density.  This  rule 
anticipates  potential  cost  savings  for  larger  versus  smaller  shelter 
complexes . 

c.  The  cnoi.ce  of  100-psi  shelter  sizes  is  100-,  500-  or  1,000-man. 

As  can  be  deduced  from  the  sample  data  shown,  a  1,000-man  shelter 
is  more  economical  than  a  one-500/f our-100-man  mix  for  sheltering 
elements  having  801  to  900  population. 
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APPENDIX  D 
SPACES  AND  COST  DATA 
FOR  MOST  EFFICIENT  SHELTER  POSTURES 
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APPENDIX  D 


SPACES  AND  COST  DATA 
FOR  MOST  EFFICIENT  SHELTER  POSTURES 

(Designed  to  Counter  a  Single  10-Mt  Surface  Burst 
on  the  Houston  At-Home  Population) 


Percent 

Minimum 

Survival 

Type 

Shelter 

Persons 

Sheltered, 

Thousands 

Spaces 

Required3, 

Thousands 

Spaces  . 
Provided0, 
Thousands 

Total 

Cost, 

Millions  of  $ 

65 

100-psi 

52 

53 

53 

20 

Fallout0 

1,174 

1,174 

1,174 

100 

Total 

1,226 

1,227 

1,227 

120 

70 

100-psi 

131 

133 

133 

51 

Fallout 

1,095 

1,095 

1,095 

93 

Total 

1,226 

1,228 

1,228 

144 

100-psi 

E| 

321 

322 

123 

80 

Fallout 

HE: 

909 

909 

77 

Total 

1,226 

1,230 

1,231 

201 

90 

100-psi 

620 

631 

635 

245 

Fallout 

606 

606 

606 

52 

Total 

1,226 

1,237 

1,241 

296 

95 

100-psi 

912 

932 

941 

367 

Fallout 

314 

314 

314 

27 

Total 

1,226 

1,246 

1,255 

393 

300- psi 

229 

232 

232 

179 

98 

100-psi 

860 

809 

902 

358 

Fallout 

137 

137 

137 

11 

Total 

1,226 

1,258 

1,271 

648 

a.  Population  in  each  grid  element  is  rounded  to  next  100. 

b.  Additional  spaces  are  provided  if  larger  shelters  are  more 
economical  for  meeting  the  requirement. 

c.  7.3-psi,  PF  =  40. 
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